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Abstract 

Fano resonance spectra for the G band in metallic carbon nanotubes are calculated as a function 
of laser excitation energy in which the origin of the resonance is given by an interference between 
the continuous electronic Raman spectra and the discrete phonon spectra. We found that the 
second-order scattering process of the q ^ electron-electron interaction is more relevant to the 
continuous spectra rather than the q = first-order process because the q = direct Coulomb 
interaction vanishes due to the symmetry of the two sublattices of a nanotube. We also show that 
the RBM spectra of metallic carbon nanotubes have an asymmetric line shape which previously 
had been overlooked. 

PACS numbers: 78.67.Ch, 73.22.-f, 42.65.Dr, 03.65.Nk 
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Raman spectroscopy of single wall carbon nanotubes (SWNTs) and graphene has pro- 
vided us with a better understanding of many optical properties which are very important 
for characterizing SWNTs and graphene not only for basic science understanding but also 
in applications Although most of the excitonic physics in the Raman spectra of SWNTs 
has been investigated intensively in terms of, for example, the excitation energy dependence 
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, chirality dependence (the Kataura plot) 



4j-|6|, Fermi energy de- 



pendence (the Kohn anomaly) |7H9|, polarization dependence 10Ml2l|. and even strain depen- 



dence 



13, 



14{, however, the asymmetric spectral shape of the G band for metallic SWNTs 



(m-SWNTs), also known as Breit-Wigner-Fano (BWF) line shape, is still not well explained 
theoretically. In a previous study, Brown et al. showed the diameter-dependent asymmet- 
ric spectral shape of the G band in which the asymmetric factor l/gswF depends on the 
density of states at the Fermi energy [lsj. Additionally, the BWF line shapes appear in 
graphite intercalation compounds (GICs) where I/^bwf depends on the staging number of 
GICs and thus also depends on the density of states at the Fermi energy 16j . Therefore, 
electrons in the gapless linear energy band of m-SWNTs should be expected to exhibit these 
asymmetry-related phenomena. 

Fano pointed out that the asymmetric feature of a broadened spectrum comes from an 
interference between a discrete excitation spectrum and a continuum spectrum jl?]]. In 
m-SWNTs, electrons in the linear energy band play an important role to give rise to the 
continuum spectra and phonons give the discrete spectra. However, the detailed mechanism 
of the BWF line shapes in m-SWNTs remains a long-standing debatable topic. Some reports 



suggest that the coupling of a col 
origin of the BWF asymmetry 



ective excitation (plasmon) with a phonon could explain the 
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20|, and some others argue that the single 



electron-hole pair and phonon coupling via the Kohn anomaly is more relevant 



partic 
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22]. 



Recently, Farhat et al. have observed a new feature of the continuum spectra exclusively 



in m-SWNTs which is ascribed to the electronic Raman scattering (ERS) 23]. The ERS 
feature (at ~500 cm -1 ) is observed in the energy region between the RBM and the G 
band and shows no phonon feature by the following arguments: (1) in comparison to the 
phonon spectral width (~l-50 cm -1 ), the ERS width is much broader (~500 cm' 1 ) and 
has a smaller peak intensity (Jers ~ 0.6/g), (2) the energy of the inelastic scattered light 
in the phonon Raman spectra is changed by changing the laser excitation energy E^, while 
the ERS peak position does not change; it keeps constant at Ma (i th Van Hove singularity 
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transition energy), (3) the ERS feature is suppressed by changing the Fermi energy, which 
indicates that the origin of this spectrum comes from electron-hole pair excitations in the 
linear band of m-SWNTs by the Coulomb interaction. 

In this Letter, we propose that the BWF feature of m-SWNTs comes from the interference 
between the G band and the ERS spectra. We calculate the exciton-exciton matrix elements 
of the Coulomb interaction which are responsible to give the ERS spectra. The exciton 
consideration is based on the fact that exciton effects in m-SWNTs are not negligible due to 
the one- dimensional carrier confinement even in the presence of the screening effect 24], |25 |. 
The calculated results of the present work suggest that the zero momentum transfer (q = 0) 
vanishes in the direct Coulomb interaction because of the symmetry of the wave-function, 
and thus a higher order Raman process is more relevant to the ERS. By considering the 
second-order Raman process, we are able to reproduce experimental results of the ERS 
spectra consistently. We will also show that the RBM spectra of m-SWNTs have a similar 
asymmetric line shape, indicating that the ERS can be coupled with both the G band and 
the RBM. 

Optical processes of the ERS consist of (i) an exciton generation via an exciton-photon 
interaction, (ii) excitation of another exciton in the linear energy band by the Coulomb 
interaction with the photo-excited exciton, and (iii) finally the photo-excited exciton goes 
back to the ground state by emitting a photon. The exciton-exciton interaction in (ii) may 
occur in a first-order or high-order process. Here, we consider up to second-order processes 
for simplicity. For the first-order process, the photo-excited exciton relaxes vertically (q = 0) 
from a virtual state \l/ vir to the Ma state after photo- absorption at a wave vector k, while 
the other exciton is created in the linear band at wave vector k' by the Coulomb interaction 
(see Fig. Ufa)). In the second-order process, on the other hand, the existence of the two 
inequivalent K and K' points in the graphene Brillouin zone leads to two different scattering 
processes, i.e. intra-valley (AV) scattering and the inter-valley (EV) scattering, shown in 
Figs. H](b) and (c), respectively. In both cases, two excitons are excited at the linear band. 
For each scattering process, we also have two cases where the two electrons at parabolic 
and linear bands may exist in the same valley (A state), or they may exist in the different 
valleys (E state). Such a symmetry labeling for A and E states is obtained from the group 
theory. After going through the electronic scattering process, the photo-excited exciton then 
returns to the ground state by emitting a photon with resonance energy E s = Ma. This is 
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the reason why the ERS peak remains at Mu even though we change the laser energy E^. 
Considering all the processes shown in Fig. [H we write the perturbed Hamiltonian as: 

H e _ e = wi±) (( k + q) > ( k ' - q) > k > k ' 

k,k',q 

X 4+q C ^_ q C k' C k' (!) 

where k and k' denote, respectively, an electron state in the parabolic and the linear band, 
while cj^Ck) is the creation (annihilation) operator in the conduction (valence) band. The 
direct (exchange) interaction K d (K x ) contributes to the two-body Coulomb interaction W 
as follows: = K d ± K x , in which + (— ) gives a singlet (triplet) state for the two 

electrons. K d and K x are expressed as j^ j: 

T^d \ ^ /~1C* (~1C* /TIC f~i\ 

A - «,k+q u «',k'-q°*,k u a ',k' 

ss'=A,B 

xKK'(q)), (2) 

tsx \ r /~ic* /~<c* /-ic 

A - «,k+q «',k'-q u «'.k «,k' 

ss'=A,B 

xft(uw(k'-k-q)), (3) 
where ^ are the tight binding coefficients for s = A, B atomic sites of the conduction 



w (a) 



(valence) band, 3?() is the real part of a complex variable, and the screened potentia 
is given by the random phase approximation (RPA): w (q) = v (q) /k(1 + v (q) II (q)) 
Here t> (q) denotes the Fourier transform of the Ohno potential, II (q) is the RPA polarization 
function, and k is the static dielectric constant due to electronic core states, a bands, and 
the surrounding material. In this calculation we used a constant k = 2.2 [6J. 

The exciton-exciton matrix element for the photo-excited exciton and another exciton in 
a linear energy band is calculated using the following formula: 

-Me t x - ex (q)=(^ / |^e-e|^ ir ) 

= ^2 ^(*k+q)c,kv^*(k'-q)c,k'v^kc,kv 
k,k' 

xH/( ± )((k + q),(k'-q),k,k'). (4) 
Here the photo-excited exciton state is defined by: 

\^ V1I ) = Y, Z ^mA C cI\9), (5) 

n,k 



(a) first-order process (b) second-order process, 



intra-valley (AV) scattering 




(c) second-order process, inter-valley (EV) scattering 




Figure 1: (Color online) (a) First-order electronic Raman process, (b) AV and (c) EV second-order 
scattering processes (q ^ 0). In both the first-order and second-order processess, the interaction 
between electrons in the parabolic band and the linear band can take place in the same valley (K 
or K') or in a different valley. 5 
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Figure 2: (Color online) Singlet state exciton-exciton matrix elements A / l c ! x _ cx calculated for a 
(15,6) tube with a diameter of 1.46 nm. Panel (a) is for the intra-valley scattering and (b) is for 
the inter-valley scattering. Label A (E) inside each panel shows the A (E) states, in which two 
electrons lie in the same (different) valley. The wave vector q is projected on the one-dimensional 
SWNT cutting lines and expressed in terms of the translational vector length T. We have T = 0.89 
nm for the (15,6) tube. 

where Z£* kv is the eigenvector of n— th exciton state solved from the Bethe-Salpeter equa- 
tion, k c and k v denote wave vectors for the electron and hole states, respectively, with 
k c = k v for a bright exciton, and \g) denotes the ground state [6j. In Eq. ([5J), we only 
use the lowest exciton state n = 0, since it gives the dominant value to the exciton-photon 
matrix element j^.The final state of Eq. (H} is given by: 



\* f ) = 5Z Z (k+q)c,kv^( k '_ q ) Cik ' v cj c c +q cj c < ;_ q c^4|( ? ). 
k,k' 
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In Fig. El we show the calculated exciton-exciton matrix elements for singlet states A4c X - cx 
in a (15,6) m-SWNT. It is noted that the matrix elements for the triplet states Aiex-cx are 
comparable with those for the singlet states, so here we only show the singlet state case. 
Surprisingly, the AV scattering matrix elements give almost a zero value at q = for both A 



and E states, as shown in Fig. [2(a). In fact, we find that the direct interaction K d vanishes 
at q = for all nanotubes. At q = 0, only the exchange interaction K x gives a small 
contribution from the AV scattering. The vanishing K d can be explained by the presence 
of three C| and one CJ coefficients in Eq. ([2]). The product of wave functions always gives 
an opposite sign when we exchange A — > B in s or s' and thus the total summation over 
A and B sublattices vanishes at q = j^. As long as we incorporate three C c and one 
C v coefficients into K d , the vanishing direct Coulomb interaction at q = is a general 
phenomenon in graphene and SWNTs systems. Furthermore, the EV scattering matrix 
elements shown in Fig. [21(b) are an even function of q because C c and C v change their signs 
by exchanging K and K' in the B sublattice, while in the A sublattice there are no changes 
in sign for C c and C v . The results from Figs. [2(a) and (b) thus imply that the first-order 
Raman process corresponding to the AV scattering at q = as shown in Fig. [2(a) makes 
only a minor contribution to the Raman spectra. Consequently, we should consider the 
second-order ERS process, in which the term in Fig. [2(a) becomes important. 

Next, to explain the Fano resonance in m-SWNTs, we calculate the Raman intensity by 
taking into account each contribution from the RBM, ERS, and G band, and considering 
all possible initial (i) and final (/) states: 

2 



2J 
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(7) 



in which we have Arbm and Aq for the phonon spectral amplitudes and Aers for the 
electronic scattering amplitude. We do not consider the C band because its position (~ 
2700 cm" 1 ) is quite far from the ERS and might not interfere with the ERS as indicated in 
Farhat's experiment (Fig. E(b)). The amplitude of each phonon spectrum can be calculated 
by: 
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(8) 



[E L - huj v -fkj B - iT v 

where v = RBM or G mode, AE m i = E^ — E m — Ei, and hu s is the scattered photon energy. 
We use a broadening factor 7 = 60 meV for the life time of the photo-excited carriers 28]. 
We also utilize the phonon spectral width for the RBM as a constant Trbm = 10 cm -1 , 
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Figure 3: (Color online) (a) Calculated result (this work) and (b) experimental results (adapted 



from Ref. 



23l |) of Raman intensity versus scattered photon energy (?iu) s ) for a (23,14) tube where 



we have the calculated = 2.10 eV and the experimental = 2.08 eV. The laser excitation 
energies E L are taken as 2.00, 2.07, 2.10, 2.14, and 2.20 eV. 

and for the G band, which consist of in-plane transverse optic (iTO) TiTo — 20 cm -1 and 
longitudinal optic (LO) IYo = 31 cm^ 1 [21]. The exciton-photon (Ai b S_^) and exciton 
phonon (Ai b)C 



cx— op 7 

,. x 13 h) matrix elements for a transition between states a — > b are taken from 



Jiang's work 



271 ] . We approximate the virtual states i — f and n = n'. On the other hand, 
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the amplitude of the second-order ERS process is given by: 
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[E^ — Tu/Ji — — hu s — 2iT x ] 



where we also consider the same virtual state approximation as in Eq. (JSJ). Here, uj% and u>2 
are the energies of the linear band excitons emitted from the exciton-exciton interaction in 
the second-order ERS process. The summation over a denotes all different processes in the 
ERS mechanism, i.e. AV and EV scattering processes. The electron-electron interaction life 
time is set at a constant value T x = 25 meV. 

In Fig. |3]^a) we show the calculated result of the E L dependence of the Raman intensity 
as a function of scattered photon energy (ftu> s ). In the present work, we only calculate the 
-El dependence of the Raman intensity near M\ 2 . Despite the two contributions from the 
LO and iTO phonon vibrations for the G mode, the splitting of the G + and G~ modes do 
not appear visibly in the spectra due to the large diameter of the (23,14) tube (d t = 2.5 
nm) which makes it has a small curvature effect. Even though we can not reproduce the 
relative intensity scale exactly from the experimental data, our calculated result can explain 
the behavior of the observed ERS as shown in Fig. E^b). The ERS feature has a very broad 
spectral width (FWHM ERS 50 meV) with a peak intensity almost comparable to that of 
the RBM. Unlike the other phonon modes, whose peak positions are shifted by changing E^, 
the ERS peak remains at the frequency of the Ma transition. At E^ = 2.07 eV, the ERS 
spectrum starts to appear and modifies the RBM and the G band line shapes. At that point, 
although El is 30 meV below Ma, the energy-momentum conservation during the exciton- 
exciton scattering may be violated by the Heisenberg uncertainty principle (At w 10 fs 
corresponding to AE w 100 meV). 

Each Raman intensity calculated from Eq. ((7j) actually gives a Lorentzian shape for all 
phonon modes and also for the ERS as presented in Fig. H^a). However, the broad feature 
of the ERS overlaps with the phonon modes and thus the interference between them gives 
rise to the asymmetric line shape, peak shifting, and the enhancement of both the RBM 
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Figure 4: (Color online) Calculated Raman spectra for a (23,14) SWNT with E L = 2.14 eV. 
The total intensity shown in panel (a) is represented by the solid line. The dashed lines show 
contributions from the RBM and G modes, while the dotted blue line is the contribution from the 
ERS. Each line shape for the RBM, the G modes, and the ERS are Lorentzian. (b) The RBM and 
(c) the G band spectra after subtracting the ERS spectrum. The dashed lines are identical with 
those in (a), (d) Fitting result of the G band spectrum to the BWF line shape of Eq. (fit)]) . Filled 
squares are calculated results and the solid line shows the BWF fitting. 

and the G bands, which can be seen in Figs. 11(b) and (c). We find that the asymmetric 
line shape of the G band after subtracting the ERS contribution clearly shows the BWF line 
shape (Fig. gfd)), fitted by 



[1 + (u - u ) fqswF^Y 



(10) 



i + [(u - u, ) AT 

where qswF, I\ and Uq are parameters to be determined. From this fitting, we can find 
and analyze the E^ dependence of the asymmetric factorl/g B wF, the spectral width T, and 
the peak position u (see Fig. [5]). According to Fano [17|, 1/<?bwf is proportional to the 
coupling constant between the continuum spectrum and the discrete spectrum. In our case, 
|1/?bwf| (FWHM or uj ) as a function of resonance condition — has a "A" ("V") 
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Figure 5: (a) Asymmetric factor (I/qbwf), and (b) spectral width and peak position of the G band 
as a function of resonance condition for the (23,14) tube. The solid and dashed arrows are given 
as a guide for the corresponding axes. . 

shape, with the maximum (minimum) peak ~ 40 meV above the resonance as depicted in 
Fig. Efc) (Fig. EDO). 1 1/9 bwf| reaches a maximum value because the intensity and the peak 
position of the ERS allows it to have a very strong overlap with the G band at that point. 
This coupling also induces the narrowing and the shifting of the G band peak closer to the 
ERS peak position. 

From our theoretical point of view, we suggest some conditions how the ERS and asym- 
metric phonon modes in m-SWNTs can be observed experimentally. Since the Coulomb 
interaction is inversely proportional to the SWNT diameter dt and also sense a curvature- 
induced band gap (~ cos8/df meV) appears for small chiral angles 9 for dt < 1 nm [29], the 
diameter range of m-SWNTs which allows us to observe the ERS and Fano resonance should 
be around 1 — 2 nm. Moreover, the finite length of m-SWNTs leads to discrete k points and 
the electron-electron interaction energy is around 60 meV; thus the nanotube length should 
be larger than 4 jum for 1 meV energy resolution. The energy of the second-order exciton- 
exciton interaction (hui) is only ~ 10 meV lower than the first-order process because the 
linear band slope is steeper than that for the parabolic band. Therefore, in order to identify 
the dominant contribution of the second-order process, the low temperature (10 — 100 K) 
gate voltage experiment must be performed. 
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In summary, we have formulated a theoretical picture of the ERS by considering the 
exciton-exciton interaction. We showed that the non-zero momentum transfer process 
gives a dominant contribution to the ERS spectra. This ERS spectrum is strongly coupled 
with the G band and the RBM and the interference with the ERS spectrum modifies the 
line shapes of the two phonon modes which results in these phonon modes having the BWF 
line shapes. The asymmetry, narrowing, and shifting of the G band induced by interference 
with the ERS are all sensitive to the peak intensity ratio and the peak distance between the 
ERS and the G band. The RBM mode also is predicted to have a similar asymmetry which 
opens up the possibility for future experimental observations and clarifications. 
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